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Using a complete basis set we have obtained an analytic expression for the matrix elements of 
the Coulomb interaction. These matrix elements are written in a closed form. We have used the 
basis set of the three-dimensional isotropic quantum harmonic oscillator in order to develop our 
calculations, which can be useful when treating interactions in localized systems. 

PACS numbers: 71.15.-m, 02.70.-c 
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INTRODUCTION 

Having an analytic expression for the Coulomb matrix 
elements is an important step for several numerical meth- 
ods, like, for example, exact diagonalization method. In 
order to describe the Coulomb interaction in three di- 
mensions, we have chosen the basis set of the isotropic 
harmonic oscillator for the single-particle wave functions, 
which, in one dimension, is written as 

V>„» = {aV^2 n *n x \y 1/2 e -5- 2 /a 2 Hnx (x/a), (1) 



where a = yj h/mu) is taken as the characteristic unit 
length. One of the reasons for the election of this partic- 
ular basis set is the Gaussian Product Theorem, which 
guarantees that the product of two Gaussian type orbitals 
(a linear combination of them in our case) centered on 
two different atoms is a finite sum of Gaussians centered 
on a point along the axis connecting them. 

In previous works, several ways to evaluate the two- 
dimensional matrix elements using different approaches 
have been studied 0, 0, @], such as restrictingto the 
lowest Landau level due to simplicity reasons [J, |a] ■ 

The purpose of this paper is to report an analytic for- 
mula for the Coulomb interaction written in closed form. 
It can be easily implemented by computer means and 
could help to improve the performance of solid state sim- 
ulations in which interactions are taken into account. 



MATRIX ELEMENTS 

In order to derive an analytical expression for the 
Coulomb interaction matrix elements we will proceed 
starting with the same approach as the one used in Ref. 
0, i.e., writing the single-electron wave function and the 
Coulomb potential as their Fourier transform integrals: 



Mr) = J MQ)e-^ dq, 

1 



V(r) = 



(2tt) 3 / 2 



V{q)e-^ r dq, 



(2) 
(3) 



where A stands for a set of quantum numbers {rii} and 
V(ri — r2) = r±2 is the Coulomb potential. Now, the 



two-particle matrix element, which, in real space is writ- 
ten as 



Kxl = J PxririMMWn-ri) 
x V>A 3 (72)^X4 ( r i) drir 2 , 
is now expressed, in momentum space, as 

;AlA2 _ 



(4) 
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i(qi - <?) 



AaAl (2tt) 
x 4L (<?2)0A 3 (<?2 + q)V(q) dq 1 dq 2 dq. 



(5) 



Eq. (jSJ) can be rewritten in a more convenient and com- 
pact form. Let us define C^,(q) and D^(q) as the fol- 
lowing convolution integrals: 
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-q) dk 


(6) 


J ip\{r)ip\ 


(r] 


e- lqr dr, 


(7) 


J ^x{k)<t>x> 


(k 


+ q) dk 


(8) 


J ^*\{r)^x 
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e lq r dr 
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Cy(-q)- 






(10) 



Substituting Eqs. JHJl and (JHJ into Eq. j£| we obtain 
1 



V 



Ai A2 
A3 A4 



(2tt) 3 / 2 



Cll (q)D^(q) V(q)dq. (11) 



Now, it is straightforward to perform the integral ap- 
pearing m Eq. 0- Using Cartesian coordinates, it is 
possible to separate all three variables and integrate in- 
dependently. For simplicity reasons, let us integrate only 
along the x variable, the result then reads: 



2<+ nl±\ 
nl\\ 2"" 



1/2 



xe 



la 2 1 A 



In — n I I 1 4 
1 x x 1 — <n — ?i 



n x — n xl /„2 2 



{a 2 q 2 x /2) , (12) 



where n\ is the quantum number referring to the i-axis 
of the particle j. We have also used the terms r£, and 



2 



n{^_, which are defined as max(n^,n* : ) and min(n^,n* : ) 



respectively. The final form for (q) will be 

c t^i^)= n c :U^- ( i3 ) 

ie{x,y,z] 

Using the relation between D and C shown in Eq. l]lUp. 
it is trivial to find out the value of the former convolution 
integral. 

It still remains to calculate the Fourier transform V(q) 
of the spherically symmetric interaction potential V(r). 



(14) 



But it will be more convenient to substitute it by 

V(q) = M (°° e-^ + <+& u du (15) 



The integration over variables q x , q y and q z can be 
performed all in the same fashion. Using the symmetry 
of the problem we only need to integrate over one vari- 
able, i.e. q x and then use the same result for q y and q z . 
Therefore, integrating over q x yields: 



u+a 2 /2)ql ^1x^ 



Wl- n i\ + Wl-nl\ 



><&t < \a\l/2)&t nl \a\l/2)dq x . (16) 



This integral does not vanish if and only if 

\ n l- n t\ + \nl = 2s x , 



(17) 



where s x = 0, 1, 2, . . .. Therefore, using the previous se- 
lection rule and the power series for the associated La- 
guerre polynomial 



k=0 



(18) 



we can write Eq. (|Po)l as 



E 



kx=0 \ x- x/ k ^ =Q x \ x- x 



(~l) fc - ( nf + 



x2 K 



2sx+2k x +2k' 



(2s x + 2k x + 2k' - 1)!! , — , x 

x (2 : + a^^/ 2 ^ (i9) 



Taking into account only the it-dependent part in 
Eq. (|19|) and its symmetric extension for y and z vari- 
ables, we end up with the last integral which will lead to 
the final result. This last integral is expressed as: 



(2u + a 2 ) 



2 \ -fi-3/2 



du 



1 1 

1 + 217^+20' 



(20) 



where Q, — s x + s y + s z + k x + k y + k z + k' x + k' y + k' z . 

Finally, collecting all the terms, we end up with the 
analytic expression for the Coulomb interaction matrix 
elements: 
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, (2s 2 - 


h 2k z + 2k' z - 


-1)!! 



22s z +fe z +fe?. 



1 + 2(s x + s y + s z + k x + k' x + k y + k'+k z + k' z ) 



(21) 
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RECURRENCE 



Eq. l|22(l can be transformed to 



Due to the six summatories appearing in Eq. 1|21|) . if 
the indices start to grow to values say, just of the order of 
tenths, the process for calculating a single matrix clement 
can be quite time-consuming, and thus, a real bottleneck 
for any numerical simulation. Using the recurrence re- 
lations that the Hermite polynomials obey, it is possible 
to find a simple iterative formula for the matrix elements 
which will accelerate the process of calculating the matrix 
elements. 

Let {n_,n + } be any pair of quantum numbers 



{it 



jk jk 



} with i £ {x,y, z} and jk £ {14,23}, satis- 



fying n + > n_. Then, the Coulomb matrix elements will 
satisfy (remaining indices ommitted for clarity) 



TT n + + 1 



V„! +1 =V n :' +2n + V„I -2n_V£_ v 



(24) 



Another interesting recurrence relation which, this time 
involves four indices {0, ti+} and {m_,m+}, is the fol- 
lowing 



— n + ,m + + l Yjn-|_+l,ro+ — n + ,m + 



(25) 



1 1 



In- - 1. 



n_-2> 



(22) 



for n_ > 0. If we consider the unnormalizcd matrix 
elements 



V. 



/r,' /r,, /r.' // ; // .' // ; 



JJ (2<<!2"^;!2<<!)^ 



iG{l,2,3,4} 
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